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Abstract 
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Blacksburg, Virginia 24061 


Abstract 


This paper is concerned with the derivation of the equations of 
motion for the Spacecraft Control Laboratory Experiment (SCULE). For 
future reference, the equations of motion of a similar structure 
orbiting the earth are also derived. The structure is assumed to 
undergo large riyid-body maneuvers and small elastic deformations. A 
perturbation approach is presented where the quantities defining the 
rigid-body maneuver are assumed to be relatively large, with the elastic 
deformations and deviations from the rigid-body maneuver being 
relatively small. The perturbation equations have the form of linear, 
non-self-adjoint equations with time-dependent coefficients. An active 
control technique can then be formulated to permit maneuvering of the 
spacecraft and simultaneously suppressing the elastic vibration. 

1. Introduction 

_Some of the contemplated NASA missions involve experiments 
consisting of the control of flexible bodies carried by the shuttle in 
an earth- orbit. Other missions involve laboratory simulations of 
similar experiments. Hence, a formulation capable of accommodating both 
types of experiments is desirable. To this end, we propose to derive 
Layrange's equations of motion for the spacecraft of Fig. 1 regarding 
the structure as orbiting about the earth and then modify these 
equations so as to describe the laboratory experiment. In the 
derivation, the shuttle is treated as a rigid body and the beam and 
antenna as flexible, distributed-parameter members. The equations of 
motion can be further modified for the case of a rigid antenna. 

The equations describing the maneuvering of a riyid space structure 
consist of nonlinear ordinary differential equations. On the other 
hand, the equations describing the small elastic displacements of a 
flexible structure relative to the riyid-body maneuver are linear 


*$upported by the NASA Research Grant NAG-1-225 sponsored by the 
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partial differential equations. Hence, the complete equations of motion 
describiny a flexible body during a maneuver represent a set of 
nonlinear hybrid differential equations. 

J Hybrid systems possess an infinite number of degrees of freedom. 

In practice, however, it is necessary to reduce the number of degrees of 
freedom to a finite one, which implies spatial discretization and 
truncation. Substructure synthesis often proves useful as a method of 
discretization and truncation, particularly in the case of distributed 
substructures. Even in the case of discrete substructures, a set of 
linearly independent vectors can be used as admissible vectors to reduce 
the number of equations of motion. 

In this paper, we propose a perturbation technique whereby the 
flexible spacecraft maneuver is assumed to consist of a combination of a 
rigid-body maneuver and small perturbations, including rigid-body 
deviations from the rigid-body maneuver and elastic vibrations. 

Regarding the rigid-body maneuver as known, the perturbation equations 
for the vibration control reduce to a set of linear ordinary 
differential equations with known time-varying coefficients. 

2 . Equations o , ’o. ion of the Spacecraft 

It is convenient to refer the motion of the spacecraft to a given 
reference frame XQygZg, where the frame can be regarded as being 

embedded in the rigid shuttle. The reference frame has six degrees of 
freedom, three rigid-body translations and three rigid-body rotations. 

We propose to derive the equations of motion by means of the 
Lagrangian approach. To this end, we must first obtain expressions for 
the^inetic energy, the potential energy and the virtual work. 

Consider.! ng Fig. 1 and denoting the position of the origin 0 of the 
frame XgygZQ by the vector R and the position of a point S in the 

shuttle relative to 0 by r, the position of S relative to the inertial 
frame XYZ. is R,. - R )• r. Moreover, denoting by a the vector from 0 to a 

nominal point A on the appendage and by u the elastic displacement 
vector of the point, the position of A in the displaced configuration 
is R n = R t a t u. It must be noted that the vectors r, a and ii an 

likely to be measured relative to axes XgygZg. In view of the above, 

the velocity of. a point S in the shuttle is 


R s = R •' «i! * r (L) 


where R Is the translational velocity and w is the angular velocity of 
the frame XqYqZq with respect to the inertial frame. Similarly, the 
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velocity of a point A in the appendage is 


R n = R + in x (a > u) + u 

vli v v ' v Ni 


(2) 


where u is the elastic velocity of the point relative to the x 0 ygZQ 
frame/ Hence, the kinetic eneryy of the spacecraft is 


'■?/ ISsI 2<llll S + 1 / iwS 
111 S m A 

- \ f |R + <d * r 1 2 clm s + \ J Hi + w x (a + u) + u| 2 dm A 

111c 111 A 


83 \ m|R| Z + \ y T l 0 “ + R* (« x S Q ) + \ f |u| 2 dm A 

ni A 

* 

+ R* [/ u dm A tuxj u dm A ] + / u*(u x a)din A 
m A n. A m A 


+ / («8 x a)- (w X u)dm A + j j |u x u | 2 dm 

m* in A 


A 


+ / H* (« * y)dm A (3) 

m A 

where 

Sq - / r dm s + / a diiu (4) 

m S n, A 

and iii^, m A and m are Liu? masses of Liu* shuttle, the appendage and the 
entire spacecraft, respectively. Also, Iy is the total mass moment of 
inertia matrix of the undeformed structure about point 0. Note that 

O 

| x | denotes the inner product x-x. 

The potential eneryy is due to the combined effects of yravity and 
strain eneryy. Assuming that the origin of the inertial coordinate 
system coincides with the center of the gravitational field, the 
gravitational potential con be expressed as 


.j 
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v t) s - Gm e [/ |R + r| _1 dm s + / |R + a + uf^imj (5) 

J m S m A 

where 11 ^ is the mass of the earth and G is the universal gravitational 
constant. 

The strain energy can be expressed as an energy inner product 
denoted by [ , ] (Ref. 1). The total potential energy then becomes 

v “ ^Ca.aJ + v s _ (6) 

The virtual work is due to external forces, including control 
forces. Denoting by f s the force vector per unit volume of the shuttle 

and by f ft the force vector per unit volume of the appendage, we can 

write the virtual work as 


« * / du s + 1 £a’ s »a d0 A < 7 > 

U S U A 

where D<j and D A are the domains of the shuttle and appendage, 
respectively. 

Next, we propose to discretize the system in space. To this end, 
we express the elastic displacements in the form of linear combinations 
of admissible functions, or 

~ y s H (a) 

where <!> is a matrix of space-dependent admissible functions and y is a 
vector of time-dependent generalized coordinates. Introducing Eg. (8) 
into Eq. (3), the kinetic eneryy takes the matrix form 


f = \ mR^K + Tj - 1 qU 3 + R^C^S^to + t; y^M^y + R^C^iy + R^C^o^i>y 

+ a 1 '* 1 ''!* + y'/ dm A y - | y r i- A 0*|)y *■ jj f L A (h.)y (y) 


where 


<I> = / 'i> dny , ■!> = / '!> a dm 


m 


'A 


m 


A 


A 


L A (w) " / I' |j) dni^, \ ( i> m 4> dm A 

m A ~ n, A 


(10a, b) 
(10c, d) 
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and 


, M ft « / * T * dm A (10e) 

' m A 

Is the mass matrix of the appendage. The symbol C represents a rotation 
matrix from the inertial frame XYZ to the x 0 yo z u frame and 1ts elements 

are nonlinear functions of Euler's angles «. The tilde over a typical 

vector v denotes a skew symmetric matrix of the form 



Recognizing that the magnitude of \i is large and u is small in 

comparison with the other vectors in Eg. (b) and ignoring terms of order 
higher than three, a binomial expansion permits us to write 


V = - Gm e [m|R| ” 1 - R* (S Q + / u dm A )|k| 3 ] (12) 

n, A 

Introducing Eq. (12) into Eq. (6) and considering Eq. (8), the potential 
energy can be written in the form 


. _ Gm m Gin 

“ V = \ * K »9 • # + , R 


% R t c t (5o + «l) 


(13) 


where 

= [*, ( 14 ) 

is the stiffness matrix of the appendage. The virtual work can be shown 
to have the expression 

6U = F T C6g + M T D(a)5cx + Q T 6g (lb) 


where 


£ - L £$ dD s + L £a dD 


A A 


A 


(16a) 
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B = L ?T £s dD S + /„ * T 4 dD A + f A* dD A3 
U S U A U A 


9 ■ /. A, d » A 


(16b) 


(16c) 


are generalized force vectors in terms of components about Xq.Yq and Zq 
and U(a) is a matrix of trigonometric functions of the Euler angles 
defined by the expression 

w = D(a)a (16a) 

Lagrange's equations of motion can be written in the symbolic form 


d / 3T\ . 3V _ Jr 


(17a) 


-A (ill . il + iii s d'm 

dt V J 3a 3a 
3a - 


(17b) 


d / 3T-, 3 T 3 V _ 

dt g(j " " 9 


(17c) 


so that, considering Eqs. (9), (13) and (15) and premultiplyiny Eq. 
(1717) by ir T , the equations of motion for the spacecraft in orbit are 


mR + C 1 S jjo + C^<I>c] + + 2C 1 1 d'y + + w 1 ) 't*g 

+ ^3 {ink + C T (S 0 + *c|) - 3[K l< T C T (S Q + +a)]J = C T F 


, T -T- 


J/.-2 , J \T_ 


IqU + w 1 1 q 0> + SyCU + 'kj + 1 [CR j-l* + J('d) h ’ii r j(fd) 1- — -^2 l 


Gin, 


Gin 


+ [w T ?> + d(w)]g + —3 ffR]S 0 = M 


Gin.. 


*• r i, 't* *T*T* Ulll 1 * 

H..y + 1> CR + ? w + / <l> m ay dm + j 1’ OR 

m i\ ~ l!i I 

v M V 

+ l a ('i>)! i + IX a (m) + L /\(‘il) + K A i.l = c i 


(Ida) 


(18b) 


(18c) 
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where 

J{aj) - j (aw + [aw])« dm A (lBd) 

in A 
and 

A . R 

B - 77 USe) 

IBI 

PremuHIplyincj Eq. (17b) by D"^ Is equivalent to considering Lagrange's 
equations in terms of quasi-coordinates (Ref. 2). Note that the 

position vector R, its time derivatives, and the Euler angles •* have 

been considered to be of arbitrary magnitude, with the result that many 
nonlinear terms appear in Eqs. (IB), 


3 • Equations of Motion for the Lab oratory Experiment 

In the laboratory experiment, the spacecraft is not actually free 
in space, but suspended from the ceiling by means of a cable or a 
beam. The following analysis applies to either case. The support is 
likely to affect the dynamic characteristics of the system. Hence, in 
the sequel, the support is added to the free model in the form of an 
elastic member. 

Considering Fig. 2, the position vector for an arbitrary point C on 
the-cable is R~ s c + w, where c is a position vector and w is the 

elastic .displacement of the cable, both of which are measured with 
respect to the inertial frame. The position vector for the point 0 is 



+ w B + e 


(19) 


where the subscript B denotes evaluation at Liu; point B and e is Use 

vector from point B (ball joint) to the point U fixed on the shuttle, 
measured with respect to the xgyQZQ frame. The velocity vector of an 

arbitrary point C on the cable is then 

R c = w (20) 

and the velocity of point 0 is 
R = w., + Hi x e 

V vlj V V 


» 


( 21 ) 
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The kinematics for the shuttle body and appendage remain the same as for 
the unrestrained spacecraft ifs space. Hence, the kinetic energy for the 
entire structure is 


; T - \ I |R C | 2 d«l c + i / |R S I V *\\ \\\ e dm A 

nu nu m rt 


m c 


'A 


3 \ I 1*1 2 dm c + \ I |w„ + a) x (e, + r) 1 2 dm s 

"'c 111 S 

\ V 

+ y J |w R + w x (e + a + u) + uj dm. 


ni,. 


3 g / lyl * 2 N lifgl + 2 ^ * 2" ^ A ^ 


m, 


1 

+ hn' (iS x § R ) + 1“ x e)* (<t> x §n) * J f |y! <Jn'r 


m„ 


+ in* / 5 din. + / u-(rn X b)dm A + / (uj x b). (w x u) dm. 
m A m A m A 

1 o 

+ tt J Is x u| dm A + / u«(s x u) dm A 


(22) 


m 


'A 


in 


A 


where 


S B 3 § 0 + me, b = a + e (23) 

The elastic displacement vectors u and w have been assumed to be small. 

I V V 


The expression for the virtual work is given by Eq. (7), but the 
potential eneryy must be modified. The acceleration of gravity will be 
assumed to be constant and can be expressed as y = - yu-, so that the 
gravitational potential is 


V 3 J (S + *H cll "r + / (!$ '• C) * U dni s >■ J (K '• 9 '• U)*y dm. 
J nu m„ 


A 


(24) 

where R is defined by Eq. (19) and is a unit vector in the l 

direction. The. elastic potential eneryy for the system can be expressed 
as 


\ 3 \ Cu,uU + \ [w,w] 


(25) 


where [w,w] is the energy inner product for the cable, which includes a 
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stiffening effect clue to the weight of the spacecraft model . Because a 
cable has little inherent bending stiffness, the stiffening effect can 
be significant, 

■' As with the appendage in the preceding section, the elastic 
displacement of the cable can be approximated by a linear combination of 
admissible functions, or 

n ( 2b ) 

where <}> is a matrix of space-dependent admissible functions and u is a 

vector of time-dependent generalized coordinates. Introducing Eqs. (8) 
and (26) into Eq. (22), the kinetic energy takes the matrix form 



+ J hiS'b dm A a + i T L A (y)y, - \ (27a) 

,n A 

where 

= / ‘I * 1 b dm A (2/b) 

in A 

has-been redefined and 

Mq = / ip r ip dnij, + (27c) 

m C 

is the combined mass matrix of the cable and of the structure lumped at 
the end of the cable, in which ^ denotes the matrix ij» evaluated at B. 

Introducing Eqs. (8) and (26) into Eqs, (24) and (2b) the potential 
energy can be written in the matrix form 


V = n *i' g + + mn + g t Eg + v, y f* A y + y 


( 28 ) 


where 


= / iji dm c 

m c 


(Z«a) 
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and 

K c = Dk^ 3 (29b) 


is the stiffness matrix of the support. Considering Eqs. (19) and (2b), 
the virtual work can be expressed as 


«w n F T C* B 6n + + Q T <$ij (30) 


where all the terms have been defined previously. The effect of the 
friction of the ball joint can be assumed in the form of an external 
torque. Hence, we let M = M + M f in Eq. (30), where M c is a vector of 

control moments and is a vector of frictional moments caused by the 

ball joint. 

Lagrange's equations remain in the symbolic form of Eqs, (1/), with 
the exception of Eq. (17a) which must be replaced by 


d f aT, , DV _ T r T, 
3t (;v) + a„ - C|i c t 

3rj 


(32) 


Using the same approach as in Sec. 2, the equations of motion for the 
laboratory experiment can be shown to have the form 


MqQ + Ip Jc T S B y + i|< B C T 0J T S B aj + ip^C T <f>y + ^JcVlJCj 
+ tjc^w 2 + io T )‘Imj + (im|»J + v ,T ) li + K c Q = 
iyS + + Sjc-^rj + <kj + [J((j) + w T J(w) 

+ [Cy]+]y * !• d(ii)]y I CCtjJSy : M 


M^c| + <l' l C'jj [3 ri + + / l j 1 id 1 Iim dm 


T-Tu 


+ , T Cg + L A (^)q + [L A (.o) t L a (.^) i K a Jij - 0 


(33a) 


(33b) 


(33c) 


where J(w) is redefined by replacing a with b in Eq. (lHd) and 
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T 'T** 

I B ■ I 0 + me 'e + e T S Q + sje (33d) 

is the mass moment of Inertia of the spacecraft about point 13, In this 
case, the set of Euler angles a is assumed to be of arbitrary magnitude, 
which is responsible for many nonlinear terms in Eqs. (33). 


4 . bimu i a Lio n _ and Con LroJ__ 

The nonlinear equations of motion for the orbiting spacecraft have 
the same basic form as for the laboratory experiment. Hence, the 
approach to the solution suggested here applies for both situations. 

Consider a first-order perturbation on the quantities It and *. 


It * R n + R. , a = a n + u n 

V ^ | w ■v [J V| 


(34a, b) 


where the first order terms and are small compared to the zero 

order terms ft Q and a^. Introducing Eqs. (34) into the nonlinear 

equations of motion, Eqs. (IB) or (33), and separating orders of 
magnitudes, we obtain zero-order and first-order perturbation equations. 
The zero-order equations can be used for the maneuvering of the 
spacecraft and the first-order equations for vibration suppression and 
rigid-body corrections. Before proceeding with this technique, we will 
first develop some expressions relatiny the perturbations in the Euler 
angles, a, with small angular deflections, (j, expressed in the body- 
fixed frame. This is done so that all the variables in the perturbation 
equations can be expressed in the body-fixed frame, the frame in which 
state measurements will be taken and actuating forces will be applied. 
Note that a set of Euler angles of arbitrary magnitude do not form a 
vector whereas the angles 3, being small, can be thought of as a vector 
quantity. 

First consider Eq. (16a) which relates the velocities of the Euler 
angles to the body fixed angular velocities, a. Introducing Eq. (34b) 

into Eq. ( 1 ba ) and neglecting higher-order terms, we obtain the 
perturbed angular velocity vector 


w 


to n + to i 

v* u ~ 1 


(3b) 


where 
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So c D ^0^0 * -1 3 D <Si)So + d ^q)Si ( 36a ,b) 

Using similar considerations, It can be shown that the body-fixed 
perturbation angles can be related to the perturbed Euler angles by the 
expression 


£ - = U(c»y)u^ {3 / ) 


Note that the vector £ Is a first-order perturbation of a set of quasi 
coordinates (Ref. 2). Taking the time derivative of Eq. (3/) and 
introducing into Eq, (36b), the perturbed angular velocity vector, w,, 
is related to the angles g by the expression 

= 10 o2 


Taking the time derivative of Eq. (3B), the perturbed angular 
acceleration becomes 


w. 



^ !• 


£ 


(39) 


Recall that the elements of the transformation matrix C from the body- 
fixed frame to the inertial frame consist of trigonometric functions of 
the Euler angles, so that a perturbation of this matrix will also 
involve the vector g. This relation can be derived using Eq. (37). 
Instead, consider the frame O' which differs from the U frame oy the 
angles (3\ Then, the transformation from the U' frame to the U frame 

is [I + 0] whore I is an identity matrix. Letting Cq be the 

transformation matrix from the 0 frame to the inertial frame, the total 
transformation matrix, C, from the O' frame to the inertial frame can be 
expressed as 


C = Cq + C x HO) 

where 

C : = 0Cq (41) 


In keeping with our objective of expressing the first-order perturbation 
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aquations in the body fixed frame, Eq. (34a) is replaced by 

B = B 0 + c [Si < 42 > 

» 

I 

where R^ is now a vector measured with respect to the 0 frame. The 

control forces and moments can also be expressed in first-order . 
perturbed form as follows: 


£ = &o + El ’ M- 3 Mq + Ml 


(43a, b) 


Introducing Eqs. (34) through (43) into Eqs. (18) and neglecting 
higher-order terms, we obtain two sets of equations of motion for the 
spacecraft in orbit. The zero-order equations, which govern the motion 
of the structure as if it were rigid, can be expressed as 


„T-T t 


Gni 


Tn / T """ U hL A aT T T 

"'Bo + c o s (fSo + c d“d‘i^o + . , 3 L'"Bo + (■ ■ ^ol'o^o-ol ° c cEo 

lio ' (M,,) 

v Qjll v 

s o c oBo + 7~73 S 0 C G^O + *0-0 + W U^Q-Q = Mq ( 44b ) 

I -o' 

The first-order equations, which govern the small scale motions of the 
structure, can be expressed as 


l"Bi '** [ w y t '*■ ^ C l0 (j *' ljj y '' ^ j ^ (jii 


.,-t l 


-2 


+ F|jO + r i>q + + fm|j i- + TT|'l'q = Fj 

* + ci^i + S W ) h “0 + ^ R -l h F n R i 

+ I qQ, + [ Iwq * h w qIq 

+ [ IqWq + oJqIqIiJq + [IqIOqIjq + S 0 H](J + i>q 

_ t v 

+ + Jjjlg + [Jq + + ili']q = Mj_ 


(45a) 


(45b) 


( J>^R 1 + E^^Ir, + <I> F [wl + ojf, + H]R, + 'T ^(5 
-1 Ol 0 0 


i 
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where 


+ [5 t s t 0 - o T 0 ]| + [? T ? 0 - j5s t 0 + ? T M]a 

V V 

+ m a 3 + L a 9 + + L A + = ^0 + 9i 


_ ^’p A A *|* I- 

HB i7T 3 [I • 3c Wo c y 

IcqI 

» * W + “3 W 

Gm 


ISio* 


= ■ [* Tc #o + ,773 *' c (£o + ?, ^o + / 4,, “fto dl ’VJ 


e rT 


•^T • 

A. * 


,T-T: 


IBo 


(45c) 


(4bd) 

(45e) 

(45f) 


For the laboratory experiment to be successful, the deflections of 
the cable (translation of the spacecraft) should be small. From purely 
rigid-body dynamic considerations, the point B (ball joint) must be 
close to the center of mass in the direction normal to the axis of 
rotation during a maneuver. Hence, the vector can be considered to 

be small in the following analysis. Because the motion of the cal. *■» ha? 
been assumed to be small, introducing Eq. (3b) into Eq. (33b), the i - 
order equations for the laboratory experiment are simply 






(46) 


The motion of the cable, ijm, can be expressed with respect to the body- 

fixed frame, as was done with the vector R-j for the orbiting 

spacecraft. Introducing Eqs. (3b), (39), (4U), (4i) and (43) into Lqs. 
(33) and neglectiny higher-order terms, the first-order equations of 
motion for the laboratory experiment can be expressed as follows: 


M c n * 2 LqM * CLq + Lq + K^]n + * ^IJ^U * 

+ b Iq + ‘i'Jt (47a) 


V I 


j 
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+ + t *B w O + “o *B + CW» 

• _ _ _ V 

4 ^B^o + w 0*B w Q + ^B®<P“0 + 


+ ?g + [wjt> + J 0 3g + [J 0 + ^J j 0 + CCq«iD*Ju 


-T 


= M n + ML + M, 

^ (J ^ U ^1 

_T M — T^T • «—T -2 

4) ijiyjj + 2«t Wq i^q + t (wq + w 0 )l'|jQ 


(47b) 



. rTq , i-rT^T 
+ T t> t L>|> u)q 

- J T 0 li * 

m 

, vT-T ,T -T 

l* ,0 u " Vo 

'■ • |,T LC 0 uJ T J|i 



V 

+ + L A g + 

V 

[U A + L A 

+ K A ] y = + 

9i 

(47c) 

where 






11 

O 

n-7 

" ’I’b^o s b-o + 

S B~0^ " 

(nW'J + /)C 0 g 


(47d) 

Mo = 

S R Cgg 




(47e) 


- [? T » 0 + I 

T-T. 

* w o b ~o 

dm A + i* T C 0 gj 


(47f) 


m ft 




For the orbiting spacecraft, Eqs, (44) can be solved for F Q ( t ) and 

Mg(t) for any desired maneuver strategy Rg(t) and Ug(t). For the 

laboratory experiment, Eqs. (46) can be solved for Mg(t) for any 

desi red ojg(t) . In either case, these quantities can then be substituted 
into the ri rst-order equations (Eqs. (45) or (4/)) producing a set of 
linear equations with known time-varying coefficients which govern the 
small deviations from the rigid-body maneuver and elastic motions of the 
structure. These linear equations can be expressed in the matrix form 


Mx + Ux + (K^ + ~ F 


where for the orbiting case 


x 

SJ 


T 




(48) 


(49a) 
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(49c,d) 
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■s/ 
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(49e) 


(49f ) 


(49cj) 


and for the laboratory experiment 


T 
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G = 
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In both cases, the mass matrix is symmetric. However, for the orbiting 
structure, the mass matrix is time invariant, whereas for the laboratory 
configuration, it is not. This is the case because the mass matrix of 
the cable, M^, must be calculated with respect to the body-fixed frame, 

and hence it is a function of the transformation matrix Cg(t). However, 

Mq and hence M can be considered time invariant if the axis of rotation 

is The axis of symmetry of the cable. In both cases, the time varying 
matrix G, multiplying the velocity vector, is skew-symmetric. This 
matrix contains the Coriolis (gyroscopic) terms. 

The stiffner^ matrix has been split into two parts, one part, K$, 

with purely symmetric terms, the other, K^, with the terms which lead 

to nonsymmetries. Of course, can be separated into symmetric and 

skew-symmetric matrices if desired. As can bo soeri i n Tc|s . (^Jy) and 
(bUf), the terms containing the angular acceleration, of the body- 

fixed frame lead to nonsymmetric terms in the stiffness matrix. In 
fact, many of these terms are skew-symmetric, arid hence, circulatory. 
Note that the perturbation of the Euler equations resulted in some 
nonsymmetric terms containing only the angular velocity vector This 

is not surprising, because it is well known that a steady-state rigid- 
body rotation about the axis of intermediate inertia is unstable. Uwing 
to the nonsymmetry of the stiffness matrices, the vibration problem is 
non-self-adjoint in both cases. 
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Because we ignored structural damping, the matrix G was found to be 
purely gyroscopic, and hence conservative. However, the stiffness 
matrix was found to contain circulatory, nonconservative terms. Hence, 
during a rotation, the structure may be vibrating about an unstable 
equilibrium state, requiring active vibration control to stabilize the 
desired maneuver. 

Gravity has a much stronger influence on the laboratory experiment, 
both because the force of gravity is greater, and because the opposing 
force is a point (cable) force rather than a body (centrifugal) force as 
with the orbiting structure. Hence, the static elastic displacement of 
the structure must be .considered in determining the stable equilibrium 
orientation. It has been previously assumed that the ball joint is in 
the vicinity of the rigid-body center of mass along the axes orthogonal 
to the cable. That is, point B is near the center of mass, or perhaps 
somewhere above it. We now wish to find the stable equilibrium 
orientation and elastic deflection of the structure for a given ball 
joint location. For the laboratory experiment, the static portions of 
liqs. (48) can be expressed as 


k stSst 


-ST 


(bl) 


where the static stiffness is 


'ST 


r ‘ ( 

0 

0 


S B^ C 0^ 


j 


<i> T C c 0 g] 


0 


(51a) 
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A J 


arid the static force is 


r -(»W'J + 'I' 1 )C Q tj 


-ST * 


s b c rfi 
_-* T C o3 


(bib) 


liquations (51) form a set of nonlinear algebraic equations which can be 
solved for x^. T , the static equilibrium orientation and elastic 

deflection of the structure. The elastic deflection is coupled with the 
angular orientation of the structure. Hence, the true center of mass 
depends on the orientation. The orientation about the longitudinal axis 
of the cable (direction of gravity) is arbitrary. This means that a 
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rotation about this axis will not affect the static orientation or 
elastic deflections. Any orientation can be made to be the stable 
equilibrium position by properly locating the ball joint using Eqs. 

(5L). Hence, a single-axis rotation, that is not affected by gravity, 
cap be accomplished about any desired axis once the proper ball joint 
location has been determined. 

The SCOLE design challenge maneuvers are time-optimal rotations 
beginning and ending in a state of rest. Provided the axis of rotation 
is the cable axis, the oigensolution of the structure before and after 
the maneuver will be the same. During the maneuver, the gyroscopic and 
stiffness matrices, and hence the eigenvalue problem, are functions of 
time. Solving the eigenvalue problem at each time step would be 
inefficient and of limited usefulness. However, a truncated set of the 
premaneuver eigenfunctions can be used as a set of admissible vectors to 
reduce the order of the equations of motion. The equations in reduced 
form can then be solved using a discrete-time technique. 

A most desirable control technique for a maneuver excites the 
elastic modes as little as possible. The obvious method is to apply a 
force that is proportional to the rigid-body mode corresponding to the 
maneuver. Because the modes are orthogonal, the other modes are not 
excited. The control forces must be distributed throughout the 
structure to maintain it in its initial state of deformation. Ul 
course, this would require distributed actuators and can only be 
approximated with discrete actuators, so that the elastic modes will be 
excited somewhat. 

For simplicity, consider a rotation about the mass center in the 
direction of a principal axis. The equation governing such a maneuver 
can~be expressed as 


M 0 (t) = a 0 (t)I 0 


(bZ) 


where a q( t ) is the desired angular acceleration, t ) is the moment 

delivering this desired performance and Ig is the mass moment of 

inertia. For reasons mentioned above, we wish to apply the moment by 
means of distributed actuators, so that Mg{t) is the resultant moment 

produced by forces F(p,t) distributed throughout the structure or 


Mg(t) = / r(p) F(p,t) dU(p) 
D 


(b3) 


Introducing the definition of Ig and Eq. (B3) into Eq. (W), the 
equation governing the maneuver can be expressed as 
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/ r(p) F(p,t) dD(p) » / a Jt) m(p) r 2 ( p ) dD(p) (54) 

D 0 U 

where r ( p) is the distance from the mass center to point p, where r(p) 
is jnormal to the axis of rotation, and m(p) is the mass density of the 
structure at point p. Considering) Eq. (54), and consistent with our 
desire for the force to be proportional to the rotational mode, the 
force density can bo expressed as 


F(p,t) = a 0 (t)m(p)r(p) (bb) 

■ 

• 

Because the force is proportional to the rotational mode, it is 
orthogonal to all of the other modes, so that it will excite only the 
desired rotational maneuver. Note that the force is also proportional 
to the mass, so that the mass acts as the control gain. This suggests 
that when this force is applied with discrete actuators, the actuators 
should be located at the point of maximum mass for maximum efficiency in 
riyid-body control. 


5. Cone! usions 

The equations of motion for the structure both in orbit about the 
earth and in the laboratory are nonlinear, even when the elastic 
deformations are small. The nonlinear terms result from the large 
rigid-body maneuver. Throuyh a perturbation approach, the nonlinear 
equations of motion can be transformed into a set of equations governing 
the riyid-body motions and a set of time-varying, linear equations 
governing small deviations from the prescribed rigid-body maneuver, as 
well, as elastic motions. The first-order equations are non-sel f- 
adje-int. The mass matrix is time invariant so that it need only be 
inverted- once in a discrete time simulation. The order of the equations 
can be reduced using the premaneuver eigenfunctions as admissible 
functions. For the laboratory experiment, a single-axis rotational 
maneuver, which is unaffected by gravity, is possible. The control 
force for the rigid-body rotation should be proportional to the 
corresponding rigid-body mode, with the inass actiny as the control 
gain. Also, the actuators should be located at the points of maximum 
mass . 
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Figure l. SCOLE Configuration 



Figure 2. Kinematical Representation of 
Laboratory Model 



